MICROCOPY  RtSOLUTION  TEST  CHART  ^ 

NAliONAi  BUR(  AU  0(  V>UNDARDS  I%.4  A 


PHOTOGRAPH  THIS  SHEET 


JO 

cx^  I 


w 

INVENTORY 


LEVEL  OK<0  Ojntyf.  f  INVhItORY 

ClmcixSctenec,  • 

E(ev«t(^  TU«e  of  Ori^fUKVcmft 

thi^wns. 

DOCUMENT  IDENTIFICATION  .  c/ay>.  7A 

lAb.  ^•.  3199-9. 


A«i^v«ttis 


DISTRIBUTION  STATEMENT 


AD  All 


ELEVATION  PLANE  ANALYSIS  OF  ON-AIRCRAFT  ANTENNAS 
Chong  L.  Yu  and  Walter  D.  Burnside 


TECHNICAL  REPORT  3188-2 
January  1972 

Contract  N62269-71-C-0296 


Naval  Air  Development  Center 
Warminster,  Pa.  18974 


ABSTRACT 


The  radiation  patterns*  In  the  elevation  plane,  of  airborne 
antennas  mounted  on  the  aircraft  fuselage  have  been  analyzed  by 
modem  diffraction  theory  Including  wedge  diffraction  theory  and 
creeping  wave  theory.  It  Is  found  that  the  fuselage  shape  has  a 
dominant  effect  on  the  elevation  patterns.  In  order  to  analyze  a 
general  convex  shape,  a  new  approach,  "section  matching  method", 
has  been  developed  In  which  the  fuselage  profile  Is  described  by  a 
set  of  points. 

This  new  method  Is  applied  to  some  general  geometrical  shapes 
which  approximate  aircraft  fuselage.  The  composite  elliptical 
cylinder  Is  the  most  ^neral  model  considered.  The  validity  of  this 
method  has  been  verified  by  modal  and  more  rigorous  6TD  solutions. 
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I.  INTRODUCTION 

It  is  well  known  that  the  performance  of  airborne  antenna 
systems  is  often  affected  by  the  aircraft  structure.  It  is  the 
distortion  of  the  antenna  patterns  by  various  structures  which  is 
the  primary  consideration  In  this  report  and  specifically  the  effect 
in  the  elevation  plane.  Based  on  crude  models  for  the  aircraft 
with  the  antenna  mounted  on  the  fuselage,  it  is  found  that  in  the 
elevation  plane  the  shape  of  the  fuselage  has  the  greatest  effect. 

Thus  the  main  concern  in  this  study  is  analyzing  the  effect  of  the 
aircraft  fuselage  on  the  radiation  patterns  in  the  elevation  plane. 

Modem  diffraction  theory  including  wedge  diffraction  theory 
and  creeping  wave  theory,  which  are  the  extensions  of  the  geometrical 
theory  of  diffraction  (6TD),  is  used  to  analyze  the  radiation  patterns 
of  aircraft  antenna  systems  in  the  elevation  plane.  This  theory 
provides  an  excellent  method  for  calculating  the  scattering  effects 
of  practical  structures  such  as  aircraft  fuselages  in  terms  of  ray 
optics.  This  approach  allows  one  to  consider  a  complicated  problem 
in  terms  of  simpler  component  parts. 

First  some  simple  basic  models  are  used  to  simulate  the  fuselage 
such  that  wedge  diffraction  techniques  can  be  applied  directly  to 
analyze  the  radiation  patterns.  In  order  to  introduce  the  curvature 
of  the  actual  fuselage,  the  theoretical  models  are  modified  to  include 
curved  surfaces.  Then  a  combination  of  the  wedge  diffraction  technique 
and  creeping  wave  theory  is  employed  to  solve  these  new  problems. 

In  order  to  include  a  more  general  geometry,  a  new  approach 
is  developed  in  which  the  fuselage  cross-section  is  described  by  a 
set  of  points  to  give  a  good  representation  of  the  actual  fuselage. 
Using  this  "section  matching  method"  the  elevation  plane  radiation 
patterns  of  general  airborne  antennas  such  as  circumferential  slot, 
axial  slot,  and  monopole  are  computed.  This  general  method  (computer 
solution)  can  be  applied  to  any  convex  body  which  can  be  described 
adequately  by  a  finite  nunber  of  points. 
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II.  THEORETICAL  BACKGROUND 

The  theoretical  technique  used  to  analyze  the  effect  of  the 
aircraft  structure  on  the  patterns  of  antennas  mounted  on  the  fuselage 
Is  basically  a  high  frequency  asynq^totlc  analysis  termed  "The 
Geometrical  Theory  of  Diffraction  (GTD)".  This  theory  was  originated 
by  Ke11er[1]  and  extensions  by  Kouyoumj1an,[23  Rudduck[3],  and  Ryan[4] 
Include  wedge  diffraction  theory  and  creeping  wave  analysis.  The 
radiation  patterns  of  antennas  mounted  on  an  aircraft  fuselage  composed 
of  sections  of  cones  and  cylinders  may  be  obtained  and  analyzed  by 
various  combinations  of  these  techniques.  Since  the  far  field  radiation 
pattern  In  the  elevation  plane  Is  of  main  Interest,  only  the  far  field 
form  of  the  electromagnetic  waves  Is  considered  In  this  study. 

A.  Wedge  Diffraction  Analysis 

The  basic  shapes  used  to  represent  the  aircraft  fuselage  are 
cones  and  cylinders  which.  In  our  two  dimensional  profile  analysis, 
reduce  to  a  wedge  type  geometry.  Thus  wedge  diffraction  analyses  can 
be  employed  to  analyze  the  radiated  fields  In  the  elevation  plane. 

The  diffraction  of  a  plane  wave  by  a  conducting  wedge  was  first  solved 
by  Sonmerfe1d.[5]  An  asymptotic  series  for  the  diffraction  function 
(Vg)  was  Introduced  by  Pauli [6]  as  a  practical  formulation  to  the 
solution  for  an  arbitrary  finite-angle  conducting  wedge.  Recently, 
Hutchins  and  Kouyoumj1an[7,8]  obtained  a  formulation  for  the  diffraction 
function  (Vg),  which  significantly  Improves  the  accuracy  over  that  of 
Pauli.  The  Improved  form  of  diffraction  function  for  a  wedge  defined 
by  [(2-n)  ir  ■  MA]  Is  given  by  (see  Fig.  1) 

0)  Vg(p,e)  -  l.,(p,e)  +  I+,(p,e) 

where  (p,e)  Is  the  field  point  and 


|r 


2 


(2) 


3 


(6) 


J-sin  (J-) 
R(e)  -  "  " 


cos  (J)  -  cos  (^) 


The  Improved  diffraction  function  has  been  prograimied  Into  a  digital 
computer  subroutine  for  nuneyical  calculations. 

The  basic  wedge  diffraction  problem  Is  Illustrated  In  Fig.  1. 


Fig.  1.  Basic  wedge  diffraction  geometry. 

With  a  line  source  located  at  (po>1*g)  with  respect  to  the  wedge  tip, 
the  far  field  form  of  the  diffracted  field  Is  given  by 

(7)  Up  (p.e)  -  \i\  [Vg  (Po,^-*o.n)  ±  Vg  (p,j,#^Q.n)] 

where  (p,e)  Is  the  far  field  observation  point  and  n  Is  related  to 
the  wedge  angle  by  HM  ■  (2-n}ir.  Note  that  the  diffraction  function 
Vg  has  been  applied  via  reciprocity  to  this  line  source  of  cylindrical 
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waves  at  the  antenna  location.  The  relationship  between  the  ray  R(e} 
and  the  field  value  U(p>e)  is  given  in  Reference  [3]  by 


(8)  U{p.e)*R(e)5 -  . 

/Zukp 

Thus,  in  terms  of  the  ray  form,  Eq.  (7)  can  be  rewritten  as 
(3)  *  Vg(p^,i|)+\j)^,n)] 


with  Vg  given  by  Eq.  (5). 

Notice  that  the  plus  sign  applies  for  the  polarization  of  the 
electric  field  perpendicular  to  the  edge  with  the  boundary  condition 


(10) 


'an 


wedge 


) 


and  the  minus  sign  applies  for  polarization  parallel  to  the  edge  with 
the  boundary  condition 


(11) 


(E 


) 

wedge 


*  0  . 


The  total  field  is  then,  using  the  superposition  principle,  the  sum 
of  the  diffracted  and  geometrical  optics  terms.  Note  that  in  Fig.  1 
cylindrical  coordinates  are  employed  with  the  z>axis  coincident  with 
the  edge  and  consequently  normal  to  the  plane  of  diffraction,  hence 
the  z-component  of  the  total  field  is  represented  by  the  scalar 
function  U.  Note  also  that  the  time  dependence  factor  (e^"^)  is 
suppressed  throughout  the  analysis. 

It  has  been  shown  in  Reference  [9]  that  the  radiation  patterns 
in  the  elevation  plane  based  on  the  two  dimensional  analysis  would 
not  be  adequate  for  an  antenna  mounted  on  finite  three-dimensional 
surfaces,  if  the  effects  of  discontinuities  such  as  wedge-type 
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junctions  and  tips  of  realistic  models  are  not  handled  with  sufficient 
accuracy.  Hence,  the  diffraction  by  curved  edges  must  be  introduced 
in  order  to  adequately  analyze  the  diffraction  effects  of  the  cone- 
cylinder  and  finite  cylinder  junctions. 

The  effects  of  edge  curvature  on  the  behavior  of  diffracted 
rays  can  be  described  using  an  extension  of  Fermat's  principle, 

"A  singly  diffracted  ray  connecting  two  points  is  a  curve  whose  length 
is  stationary  among  all  curves  connecting  these  two  points  and  having 
one  point  on  the  edge."[l]  By  applying  conservation  of  energy  within 
the  astigmatic  energy  flux  tube  of  diffracted  rays  depicted  in  Fig.  2, 


Fig.  2.  An  astigmatic  ray  tube. 


one  obtains  the  relation 


(12)  ^o^**''o  “ 


where  da^  is  the  cross  section  of  the  flux  tube  at  reference  point 
0,  having  principal  radii  of  curvature  p-j  and  p^;  and  do  is  the  cross 
section  at  t,  having  principal  radii  of  curvature  p<|  +  £  and  P2  £• 
A^  and  A  are  the  field  amplitudes  of  the  diffracted  rays  at  point  0 
and  i.  It  follows  from  elementary  geometry  that 


(13) 


do©  (p^+£)(p2+^) 
do  P|P2 
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Inserting  this  result  Into  Eq.  (12),  yields 


04) 


A  .  I  Pl^2 _ 


or 


This  Is  the  equation  which  governs  the  behavior  of  all  rays  (Incident, 
reflected  or  diffracted).  The  expression 


Is  the  spatial  attenuation  factor  or  "ray  divergence  factor,  (RDF)". 

It  Is  Interesting  to  note  that,  when  t  «  -p^  or  -P2»  the  field 
amplitude  becomes  Infinite  and  the  ray  optics  field  Is  Invalid.  The 
congruence  of  rays  at  the  lines  1-2,  3-4  Is  called  a  caustic.  It  has 
been  demonstrated  that  In  the  region  away  from  the  caustics  the  GTD 
solution  Is  correct  If  a  phase  shift  of  (-ir/Z)  Is  Introduced  upon 
traversal  In  the  direction  of  propagation  through  a  caustic  line. 

Now  consider  the  rays  from  a  line  source  Incident  on  a  curved 
surface  as  shown  In  Fig.  3.  It  Is  known  that  If  the  Incident  rays 
strike  the  edge  obliquely,  making  an  angle  0  with  the  curved  edge, 
the  diffracted  rays  lie  on  the  surface  of  a  cone  whose  half  angle 
Is  equal  to  0.[1O]  It  Is  clear  that  one  of  the  caustics  Is  at  the 
edge  and  the  other  one  Is  located  at  a  distance  p  from  the  edge.  The 
distance  p  between  the  edge  caustic  and  the  second  caustic  can  be 
found  according  to  Kouyoumj1an[11]  by 


(16) 


1 

p 


A  A  A 
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Fig.  3.  Vector  relations  for  diffraction  by 
a  curved  surface. 


xhere 


Pq  ■  the  radius  of  curvature  of  the  edge  at  Q 
1,  d  ■  directiwis  of  Incidence  and  diffraction 
n  >  the  normal  vector  of  the  edge  . 


The  angle  0  Is  easily  found  as 


(17)  cos  6  *  i  •  t  0  £  e  £  ir 

where  t  Is  the  tangent  vector  of  the  edge  with  the  tangent  and  normal 
vectors  In  the  plane  of  the  figure.  Note  that  the  angles  e*  and  6 
are  given  by 

(18)  -sin  6'  =  1  •  n  -ii/2  £  e'  £  3Tr/2 

A  A 

(19)  COS  6»d*n  -ir£6£ir. 

Equation  (16)  In  combination  with  rqy  optics  allows  the  determination 
of  the  effects  of  edge  curvature  upon  the  diffracted  rays. 

Figure  4  Illustrates  the  diffraction  of  a  "tube"  of  roys 


Fig.  4.  Diffraction  by  the  end  of  a  finite 
circular  cylinder. 

originating  on  the  surface  of  the  cylinder  at  the  antenna  location 
and  which  propagate  to  the  end  of  the  cylinder.  The  caustic  distance 
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of  these  diffracted  rays  is  given  by  Eq.  (16).  Applying  Eq.  (16)  we 
obtain 


(20) 


11.  sine 
02^2  *  * 


where  p2  is  the  caustic  distance,  as  illustrated  in  Fig.  4.  One 

notices  that  the  above  equation  predicts  that  the  caustic  distance  is 

infinite  for  values  of  e  =  =  sin“^(~)  such  that 

c  x>2 


(21) 


1_ 


a  +  sin  6, 
_ £ _ c 

a 


=  0  . 


A  similar  expression  can  be  derived  for  planes  removed  from 
the  s  0  plane,  in  which  case  the  caustic  distance  is  dependent 
upon  both  e  and  However,  only  the  (^  -  0  plane  is  of  interest  in 
this  study.  Note  that  for  d  <  the  caustic  distance  becomes 
negative  in  sign.  This  caustic  is  shown  in  Fig.  5.  The  significance 
of  the  negative  value  of  P2  is  that  the  diffracted  rays  pass  through 
the  caustic  location  before  proceeding  to  the  far  field. 

The  diffracted  field  from(Dcan  be  calculated  by  using  ray 
optics,  conservation  of  energy,  and  wedge  diffraction  theory.  The 
form  of  the  three  dimensional  field  incident  upon  the  edge  is  given 
approximately  by 


“jlct2  “jlci>2 

(22)  hJ  »  F(0-O.  ♦•0)  F(e-0.  ♦•0)  ^  ^ 

^  c  *^^2 

where  F(e,4)  is  the  far  field  pattern  of  the  antenna  in  the  direction 
(e,4)  (see  Appendix  V  in  inference  [9].)  Note  that  only  the  case  of 
the  magnetic  field  tangent  to  the  perfectly  conducting  cylinder  is 
considered  here  since  the  boundary  condition  on  the  perfectly  con¬ 
ducting  surface  requires  that  the  tangential  electric  field  vanishes 
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Fig.  5.  Caustic  caused  by  the  Intersection 
of  the  diffracted  rays. 


on  the  surface,  precluding  the  propagation  of  a  tangential  electric 
field.  The  factor 

-ikty 
e _ 

outside  the  brackets  In  Eq.  (22)  corresponds  to  the  form  of  the 
Incident  cylindrical  wave  field  for  which  the  two  dimensional  dif¬ 
fraction  function  Vg(i2»  ’’g)  accounts  for  the  diffraction  by 
the  discontinuity  In  the  ^*0  plane.  This 

-jki- 
e _ 
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field  dependence  Is  contained  In  the  functional  expression  for 
VgU2t  ir-e>  02)  and  thus  this  factor  appears  1np11c1t1y  when  Vg  Is 
used.  The  diffracted  magnetic  field  H^2^,  Is  now  expressed  using  the 
wedge  diffraction  coefficient  and  ray  optics  as 


(23) 


.0) 


HJ''(8.^*0)  =  [H^  F(e.O.  ♦=0)/y^  •  •  Vg(t2.  ir 


P2 


-jkr. 


ng)  •  - — 

^  /r~ 

/rg 


The  ray  divergence  factor 


1  P2 

\| 

accounts  for  the  divergence  (or  convergence)  of  the  diffracted  rays 
In  the  plane  which  contains  the  ray  path  and  Is  orthogonal  to  the 
<^■0  plane.  Inserting  the  expression  pg  Into  this  square  root  gives 


(24) 


RDF., 


atn 


aip^^Ta+ipTneT 


In  the  far  field  (rg-^*)  It  becomes 


(25) 


RDF. 


2  ^ 


a£« 


rglangSlne) 


It  Is  apparent  that  for  Eq.  (25)  Is  singular.  Radiation  at  the 
angle  e>e^  at  which  the  ray  optics  picture  falls,  produces  a  caustic 
of  the  diffracted  rays  in  the  far  field.  This  caustic  has  the  fom 
of  a  parallel  "bundle"  of  diffracted  rays.  The  diffracted  field  In 
the  caustic  region  can  be  analyzed  by  the  equivalent  diffraction 
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currents  technIqueCA]  which  is  outside  the  scope  of  this  study.  Thus, 
only  the  ray  divergence  factor  In  the  Illumination  region,  I.e.  O^e^ir, 

Is  encountered  In  this  analysis  and  that  In  the  shadow  region  Is  approxi¬ 
mated  by  the  previous  expression  with  e  set  equal  to  ir.  This  approxi¬ 
mation  gives  good  results  In' the  prediction  of  radiation  patterns  In  the 
elevation  plane  as  will  be  shown  later. 

B.  Creeping  Wave  Analysis 

When  an  Incident  ray  strikes  a  smooth,  curved  perfectly 
conducting  surface  at  grazing  Incidence,  I.e.,  at  the  shadow  boundary, 
a  part  of  Its  energy  Is  diffracted  Into  the  shadow  region.  To 
describe  this  phenomenon  Ke11er[12]  Introduced  a  class  of  diffracted 
rays  which  Is  now  well  known  as  creeping  waves.  These  r«y  paths 
Include  the  points  and  Q2  which  form  a  curve  on  the  diffracting 
surface.  However,  the  actual  concept  of  creeping  waves  was  Introduced 
by  Franz  and  0eppennan.[13,14]  A  basic  creeping  wave  concept  Is 
presented  In  the  following  discussion.  The  context  presented  below 
Is  basically  taken  from  "Asymptotic  High-frequency  Methods"  by 
KouyoumjIan.CiO] 

The  diffraction  by  a  smooth  curved  surface  Is  shown  In  Fig.  6 
In  which  0  Is  the  source  point  and  P  Is  the  observation  point  In  the 
shadow  region.  Applying  Fermat's  principle,  the  line  0Q^Q2P  Is  the 
shortest  distance  between  0  and  P  which  does  not  penetrate  the  surface. 

In  detail,  a  ray  Incident  on  the  shadow  boundary  et  Q-j  divides;  one 
part  of  the  Incident  energy  continues  straight  on  as  predicted  by 
geometrical  optics,  a  second  part  follows  the  surface  s  Into  the 
shadow  region  as  a  surface  ray  shedding  diffracted  rays  tangentially 
as  It  propagates.  With  t,  n  and  b  being  the  unit  vectors  in  the 
direction  of  Incidence,  normal  to  the  surface  s  and  blnonaal  to  the 
surface  (b  ■  t  x  n),  respectively,  the  Incident  field  E^(Q^)  may  be 
resolved  into  Its  normal  and  tangential  components  (n  •  C^(Q^}  and 
b  *  ^(Q^)).  It  Is  assumed  that  these  two  components  Induce  surface 
ray  fields  which  propagate  Independently  of  each  other  along  the 
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Fig.  6.  Diffraction  by  a  smooth  curved  surface. 

geodesic  arc  between  and  Q2.  From  Reference  [10]  the  surface  r^y 
field  Ae^*^  at  is  related  to  the  incident  fields  at  by 


A(Qi)  e‘ 


Ds(Qi)  V  rCQ,) 


where  D5(Q’|)  is  the  scalar  diffraction  coefficient  for  a  soft  surface. 
The  amplitude  of  the  surface  ray  is  assuned  to  be  governed  by  the 
conservation  of  energy  between  a  pair  of  adjacent  surface  rays.  Hence, 
the  amplitude  behavior  of  tiie  fields  is  given  as 

^2 

_  [-/  «(f)df] 

dni  Qn 

(27)  * 


where 


dn^  and  dn^  *  the  separation  between  a  pair  of  rays  at 
Ql  and  Q2>  respectively. 

a(t)  «  the  attenuation  constant  which  Is  a  function 
of  t,  the  coordinate  along  the  surface  ray, 
because  1t  depends  on  the  local  radius  of 
curvature . 

The  attenuation  constant  a  Is  Introduced  due  to  the  tangential  shedding 
of  rays  as  the  surface  ray  propagates.  It  Is  seen  from  Fig.  6  that 
Q2  Is  a  caustic  of  the  diffracted  field  and  the  second  caustic  Is 
located  at  a  distance  p  from  Q2.  Thus,  the  tangential  component  of 
the  diffracted  field  which  radiates  from  Q2  towards  P  can  be  found, 
as  In  the  previous  edge  diffraction  case,  with  one  of  the  caustics 
used  as  a  reference  point  and  Is  given  by 


(28)  b2  •  «  05(02)  A(Q2) 

From  Eqs.  (26),  (27)  and  (28)  there  results 


It  Is  found  that  b^  *  I^(Q^)  excites  an  infinity  of  surface 
ray  modes  each  with  Its  own  diffraction  coefficient  and  attenuation 
constant.  Thus,  the  expression  In  Eq.  (29)  Is  replaced  by 
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(») 


i,  ■  fcp)  -  CV  .-JClc(t«)] 

I2 


I  “sJOl)  “smOz)  ' 

m 


■I  a,(f)df 


Equation  (30)  relates  the  diffracted  field  at  P  to  the  incident  field 
at  for  the  soft  surface  boundary  condition. 

An  expression  similar  to  Eq.  (30)  is  also  obtained  for  the 
normal  component  of  the  incident  field;  in  this  case,  the  scalar 
diffraction  coefficients  and  attenuation  constants  for  the  hard  surface 
replace  those  of  the  soft  surface.  Therefore  the  vector  diffracted 
field  at  P  can  be  written  in  terms  of  the  electromagnetic  field 
incident  at  as 

(31)  r'(P)  -  Cn'j  n",  v(1.2)  +  bj  b,  U0,2)]  •  E^O) 

in  which  v(l,2),  u(l,2)  are  equal  to 


dni 


.-jkt 


m 


V) 


m 


(f)  df 


with  the  subscripts  h,  s,  respectively,  added  to  and  Note 
that  and  Q2  have  been  replaced  by  1  and  2  for  the  sake  of  brevity. 
Finding  dn-j ,  dn2.  and  p  is  simply  a  matter  of  differential  geometry 
involving  the  rays  and  surface;  this  is  discussed  at  length  in  Levy 
and  Keller.[12]  The  generalized  diffraction  coefficient  and  attenu¬ 
ation  constant  can  be  found  in  Reference  [15]. 

The  diffraction  thus  far  discussed  is  applied  to  the  open 
curved  surface.  For  a  closed  surface,  each  surface  ray  mode  produced 
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infinite  nuirber  of 
n-th  encirclement 
sed  surface.  Thes 


region  Is  the  sun  of  these  two  diffracted  fields  from  and  Q^.  The 
total  field  at  any  point  In  the  Illuminated  region  Is,  by  the  super¬ 
position  principle,  the  sum  of  Incident,  reflected  and  diffracted 
field.  A  detailed  discussion  of  this  subject  can  be  fbund  In 
Reference  [15]. 

III.  ELEVATION  PATTERNS  FOR  SIMPLE  FUSELAGE  MODELS 

The  fuselage  Is  first  approximated  by  some  simple  models 
composed  of  cylinders,  spheres  and  cones  as  shown  In  Fig.  8.  These 
models  are  simple  approximations  to  the  fuselage  In  which  the  vertical 
stabilizer  Is  Ignored.  More  general  approximations  and  the  effect  of 
the  curvature  of  aircraft  fuselage  will  be  considered  In  a  later 
section. 

Ryan[9]  has  demonstrated  that  the  radiation  pattern  of  an 
antenna  at  an  arbitrary  location  on  a  finite  circular  cylinder  can  be 
treated  as  a  two-dimensional  problem  through  the  use  of  three  dimensional 
ray  optics,  wedge  diffraction  and  creeping  wave  techniques.  It  also 
has  been  shown  that  the  diffraction  by  a  three  dimensional  cone  tip 
can  be  approximated  as  If  the  tip  were  a  two-dimensional  wedge. 
Especially  In  the  elevation  plane,  this  approximation  yields  fairly 
good  results.  Thus,  In  the  following  development  we  shall  use  similar 
approximations  to  analyze  the  radiation  patterns  In  the  elevation  plane 
for  an  antenna  mounted  on  the  surface  of  an  aircraft  fuselage. 

In  order  to  Illustrate  the  application  of  the  wedge  diffraction 
analysis  let  us  consider  the  model  In  Fig.  9  where  a  circumferential 
slot  Is  mounted  on  the  surface  at  a  distance  d^  from  edge  1  and  (a-d^) 
from  edge  5.  This  model  has  a  cone  capping  on  a  finite  circular 
cylinder.  Thus,  the  hard  surface  boundary  condition  Is  employed.  The 
majpietic  field  radiated  by  the  antenna  Is  grazingly  Incident  on  edge  1 
and  edge  5  causing  single  diffractions  at  these  edges,  which  In  turn 
gives  rise  to  the  double  diffractions  at  edge  2  and  4.  Subsequently 
triple  diffraction  occurs  at  edge  3  due  to  the  doubly  diffracted  field 


Fig.  8.  SliKple  models  for  fuselage. 


(a) GEOMETRY  OF  MODEL 


(b) RAY  geometry 


(C)  PHASE  CORRECTION 
Fig.  9.  Simple  geometry  of  an  aircraft  fuselage. 


Incident  on  the  discontinuities.  The  Interaction  between  edge  3  and 
4  produces  higher-order  diffracted  terms.  However,  the  effect  Is 
so  small  that  one  can  neglect  the  higher-order  terms  without  causing 
any  significant  change  In  the  radiation  patterns.  In  terms  of  the 
far  field  pattern  F(e,«>)  of  the  source  the  radiated  field  In  the 
elevation  plane  Is  given  by 


-jkrs 

(32)  H.(0.*»O)  -  [H  F(e.*=0)//F  ]  ^ 

where  0  ±  0  i  180® 

H  -  a  normalization  factor  which  Is  the  magnetic  field 
0 

Intensity  at  the  pattern  maximum  at  r^  *  1  unit,  as 
projected  from  the  far  field. 

Note  that  the  factor  e^*^^  Is  assumed  understood  and  Is  suppressed 
throughout  this  report.  From  the  relationship  given  In  Reference  r3], 
which  relates  the  field  value  H(p,e)  to  the  ray  R(e),  one  cen  express 
the  source  field  In  ray  form  as 


(33) 


H5(r,0) 


«0 

g-jkr  rRs(8) 

/r 

/r  ’  L  ^  . 

-J(kr+J) 
e _ 


0  <  0  <  180® 


where  R  (e)  Is  the  Incident  ray  and  Is  given  as 


(34) 


85(0)  ■  [Hq  F(e)]  /RR"  e 


0  <  e  <  180* 


Now  If  we  normalize  this  equation  by  the  quantity 


i 
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then  the  Incident  ray  Is  simply  expressed  by 

(35)  RjCe)  »  F{e)  Of.eiiso*’  . 

This  simplifies  our  analysis,  consequently  the  ray  form  of  the  radiation 
field  Is  employed  throughout  the  following  development. 

The  singly  diffracted  rays  In  any  direction  resulting  from 
diffractions  by  the  Incident  rays  striking  the  discontinuities  at 
edges  1  and  5  are  given  at  edge  1  by 

(36)  R{l^(e)  =  Vg  (d^.  180“  -  6,  n^)  •  RDF^ 

0“  <  e  <  180“ 

r{^^(0)  =  [Rs(e*0“)//?^]  Vg  (d^.  540“  -  6.  n,)  •  RDF^ 

360“  -  1  0  <  360“ 

and  at  edge  5  by 

(37)  R^^^e)  =  [Rj.(0»18O“)//?:H5]  Vg  (a-dg,  e,  ng)  •  RDFg 

0“  <  e  <  360“  -  og 

where  >  the  wedge  angle  of  edge  1 

ag  *  the  wedge  angle  of  edge  5 
np  ng  -  the  constants  related  to  the  wedge  angles  a-j 
and  og,  respectively,  by  the  relation 
HA  ■  a  ■  (2-n)Tr 
RDF  ■  ray  divergence  factor. 

The  expression  for  the  RDF  can  be  easily  obtained;  however, 
it  has  been  modified  to  avoid  the  singularity  Involvement  due  to 
the  caustic  effect  as  discussed  In  Section  I  and  Is  given  by 


0<e<i80* 

iso^ieiaeo* 


(39)  RDFg  * 


This  modification  has  been  proved  satisfactorily  In  the  radiation 
pattern  prediction  as  one  shall  see  later.  Note  that  the  subscript 
denotes  the  edge  and  the  superscript  denotes  the  order  of  the 
diffraction. 

The  doubly  diffracted  rays  at  edges  2  and  4  are  produced 
when  the  singly  diffracted  rays  and  Rg^  from  edges  1  and  5, 

respectively,  arrive  at  edges  2  and  4.  In  addition,  the  singly  dif¬ 
fracted  rays  from  edges  1  and  5  travel  In  the  opposite  direction  and 
strike  edges  5  and  1,  respectively,  causing  another  set  of  doubly 
diffracted  fields.  These  doubly  diffracted  rays  can  be  obtained  In 
a  similar  manner  and  are  given  at  edge  2  by 

0<(9£l80*-B2 


lW+B2<e<360* 
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(40) 


R{’^(e»360“-B2) 


R{^^(e-360*-B2) 


e(a-d5) 


Ue+(a-dj)s1ne]  r 


a-d. 


S  r 


0<6<180'’ 


180“<e<360*. 


at  edge  4  by 


(41) 


r|^^(6)  = 


R|^^(e=360“-a5) 


'^B^*^4*®~®4**™4^  *  *^^^4 


a^e<360“-o^ 


at  edge  1  by 


(42) 


Rf^(0)  VB(a.l80“-e.n,)  •  RDFg, 

0<e<180" 


®  ♦  Vg(a  ,54O®“0  ifi'^ )  • 


RDF^ 


18O*+a^<0<36O® 


and  at  edge  5  by 


(43) 


^(2)< 


Rr^(e) 


~  •  Vg(a,0,n5)  •  RDF^j 


0<e<360®-ae 


where 

(44) 

Rjl  »  R^’^(e-0*) 

(45) 

R^5  -  Rj^^(e«180“) 

(46) 

a  . 

®2  T 
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(47) 


RDF, 5  -  ROFfil 


e  a 

(e-<-a  sine)  r 


IF 


0<e<l80» 

180*<e<360“ 


Notice  that  the  ray  divergence  factor  RDFg,  at  edge  1  due  to  the 
arriving  of  ray  Rg,  is  different  from  RDF,  sinply  because  of  the 
distance  it  travels.  For  the  same  reason,  RDF,5  is  different  from 
RDFg.  However,  RDF,^  is  equal  to  RDFg,  because  of  the  symmetry  of 
the  problem.  Since  the  dimensions  of  the  fuselage  are  large  compared 
with  the  wavelength  the  effect  of  the  RDF  at  edge  4  is  thus  insigni¬ 
ficant  and  can  be  neglected,  i.e.,  o*ir  or  RDF^  obvious 

that  t^e  ray  divergence  factor  at  edge  2,  RDF,,  is  essentially  equal 

topf- 

Similarly  the  triply  diffracted  ray  at  edge  3  can  be  found  and 
written  as 


(AS) 


O£0<P2 


R|^^(e-180“+B2) 


Vg(d23,360*+B2-e,n3)  •  RDF^ 
18O"<0<36O“ 


As  we  mentioned  earlier,  the  dimensions  of  the  aircraft  fuselage  are 
large  and  also  the  contribution  of  the  third-order  diffractions  in  the 
far  field  radiation  pattern  are  very  small,  in  which  case  the  effect 
of  the  RDF  at-ieilge  3  is  negligible.  Thus  0  in  RDF^  is  again  set  to 
T  or  ROFj  ■  .  Here  we  have  also  neglected  those  contributions 

of  the  higher-order  terms  which  are  due  to  the  multi -interactions 
among  the  edges,  since  they  have  little  effect  in  the  radiation  pattern. 
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with  the  phase  reference  point  at  the  point  Q,  the  mid-point 
of  the  fuselage  as  shown  In  Fig.  9,  the  rays  radiated  from  the  source 
and  the  edges  are  multiplied  by  appropriate  phase  correction  terms, 
giving 

-jkq  cos(6.+e) 

(49)  R(0)  =  i^e  ^  ^  0<e<180“ 

*  /r 

(50)  R,(9)  .  [5(9=0)  •  VB((lj,180»-9.n,)  +  F(9.180»)  ■ 


Vg(a,180®-e,n^)  ,  (e+as^nejr 


Jkesine 


0<e<180‘* 


=  F(e=0)  •  Vg(d 


s* 


540“-e,n^)  +  F(e=180“)Vg(a-d5,0‘’,n5) 


•  Vg(a.540®-e.n^) 


g jkesine 


180“<e£360'» 


(51)  R2(0)  »  [f(0-O»J  •  Vg(d5,180“+B2,n,)  Vg(d,2*’®°‘’-V®»"2^  ’ 

gjkbcose 

0<e£l80‘’-B2 


•  F(0-O*)  Vg(dg,180»+B2.n,)  •  Vg(d,2»®^“®2'®»''2^ 


gjkbcose 
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(52) 


R3(0)  »  F(6=0‘‘)  Vg{d5,180“+B2.n,)  •  VgCd^g.SeO^-ZSg.ng)  * 


(53) 


(54) 


where 


-jkesine 


O<0<p2 


F(0=O»)  Vg(dg,180‘‘+62.n,)  • 


'^B^^23*^®°'*’®2"®’"3^  •  — 


-jkesine 


180“<e<360“ 


R4(0)  *  F(0-18O“)  •  Vg(a-dg.360*-o5,ng)  •  Vg(dg^,e-a^,n^)  • 

-jkq^cos(6-04) 


e_ 

/r 


Oii<0<36O®-O/| 


1)5(6)  =  [F(e.l80")  •  Vj(,.d^.6.n5)|^;^5^;;Z* 


F(6.0*)  •  V5(d,.0',n,)  .  V5(..6,n5) 
jkq5COs(0+B5) 


§_ 

✓F 


0“<e<180* 


R5(0)  -  [F(o-180'»)  •  Vg(a-d3.0.n5)  +  F(0»O‘»)  •  Vg(d^.O“,n,) 


VB(a,e.n5)]  •  — 

/r 


-jkq5Cos(e+B5) 


180*<e<360*-c 
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(55) 


Pg  =  tan‘^(|) 

»  \ld^+e^ 

q4  -  qg  *  JaV 

84  *  85  »  tan-’{|) 

a^.ag  -  the  wedge  angle  of  edge  4  and  5,  in  our 
case  here  o^*ag*90®. 

The  total  field  in  ray  form  at  any  observation  point  then  is 
simply  the  superposition  of  the  direct  ray  and  diffracted  r^ys.  Hence 
the  radiation  pattern  in  the  elevation  plane  is  obtained  as 


(56) 

8j(0)  ~  Rj (q)'*^^ (0 )^82(® ^ 

0^6<32 

(57) 

8^(6)  =  R5(0)+R,(e)+R2(e)+R5(e) 

B2<0<9O* 

(58) 

Ry(0)  *  Rj(0)^R^ (0)^R2(0)'^R4(0)'^R5(0) 

90*<e<180*-B2 

(59) 

Rt(8)  •  Rs(0)+Ri(8)+R4(0)+R5(0) 

18O"-B2<0<18O* 

(60) 

Rt(0)  -  R3(e)+R4(0)+R5(e) 

180“<e<180“+B2 

(61) 

8^(0)  »  82(0 >+83(0 >+84(0 >+85(0) 

18O*+62i0l27O* 

(62) 

8^(0)  -  R2(0)-*-R3(0)+R4(0) 

270“<e<360*-82 

(63) 

Rj(0)  ■  R^(0)+R2(0)+R3(0)+R4(0) 

36O*-B2<0<36O* 

A  calculated  elevation  pattern  for  the  model  of  Fig.  9  with 
a  circumferential  slot  located  at  d>8''  from  the  edge  1  is  shown  in 
Fig.  10,  along  with  a  measured  pattern. [9]  It  is  seen  that  the 
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f*e.95e  6Hf 


Fig.  10.  Elevation  plane  patterns  with  a  circumferential 
slot  mounted  on  a  simple  fuselage  model. 

theoretically  calculated  radiation  pattern  compares  very  favorably 
with  the  measured  result.  It  Is  noticed  that  the  field  pattern  has 
more  coverage  In  the  forward  region  than  that  In  the  backward  region, 
this  Is  simply  because  the  source  Is  located  In  the  front  portion  of 
the  model.  As  discussed  In  Reference  [16],  the  pattern  has  more 
ripples  In  the  fbrward  direction  and  less  In  the  backward  direction, 
which  Is  also  due  to  the  location  of  the  antenna  system  with  respect 
to  the  diffracting  edge. 


The  model  shown  in  Fig.  11  is  a  fication  to  the  model  des¬ 
cribed  above*  in  which  another  cone  is  mounted  on  the  other  end  of  the 
circular  cylinder.  In  general*  the  sizes  of  these  two  cones  will  be 
different:  however*  In  the  following  development  it  Is  assumed  that 
they  are  dimensionally  the  same.  Due  to  the  geometrical  similarity, 
the  ray  components  necessary  to  describe  the  radiation  pattern  of  this 
model  are  similar  to  those  in  Eqs.  (56)  through  (63)  and  are  easily 
obtained  with  only  some  minor  modifications.  For  example,  the 
expression  for  the  doubly  diffracted  field  in  ray  form  from  edge  5  is 
given  as 


-jk(a+c)cose 

(64)  R^‘^^(e)  *  R^'^(18O“+05)  •  Vg(dg5*e-05*n5)  •  ^ 

*  F(e«180*)  •  Vg(a-d,l 80^+35 ,ng)  ’  ' 


where 


--jl((a+c)cose 

✓r 


35<6060*-35 


H 


the  wedge  angle  at  edge 


5. 


As  in  the  previous  case*  we  have  made  some  approximations  which 
have  simplified  our  analysis  without  causing  any  significant  error. 
Hence,  the  radiation  pattern  is  obtained  by  simming  all  of  the  rays 
which  contribute  in  the  various  regions  and  is  given  by 


(65)  Ry(6)  ■  Rj{e)+R.j (e)+R2(8)+Rj(6)+Rg(e) 

(66)  Rj(9)  ■  Rj(8)+R>j(e)+R2(o)+R5(o)’*'R5{6)  o£e£l80*-o 

RT(e)  •  R5(9)+R,(e)+R4(e)+R5(e)-»-R5(e) 


(67) 


180*-o<e<180 


(o)  GEOMETRY  OF  MODEL 


(c)  PHASE  CORRECTION 


Fig.  n.  Simple  geometry  of  aircraft  fuselage. 


(68) 

RT(e)  = 

R3(0)+R4(0)+R5(0)+R5(e) 

180“<e<180“+a 

(69) 

RT(e)  = 

R2(o)+R3(e)+R4(e)+R5(0) 

18O“+o<0<36O“-a 

(70) 

RT(e)  » 

Rl(0)+R2(e)+R3(0)+R4(0) 

360‘‘-a<e<360“ 

where 

a  =  og  =  ttg  =  wedge  angle  of  edge  2  or  5. 

The  calculated  radiation  pattern  in  the  elevation  plane  for  a 
circumferential  slot  mounted  on  the  model  of  Fig.  11  at  a  distance 
d>8"  from  edge  1  is  presented  in  Fig.  12.  Compared  with  Fig.  10,  one 
sees  that  the  backlobes  in  the  deep  shadow  region  are  lower  for  this 
model.  The  reason  is  apparent,  since  the  addition  of  the  extra  cone 
has  reduced  the  scattering  effects  at  edges  4  and  5  and  thus  reduced 
the  level  of  backlobes.  The  calculated  radiation  pattern  in  the 
elevation  plane  obtained  by  replacing  the  circumferential  slot  by  a 
1/4x  monopole  is  shown  in  Fig.  13.  In  order  to  see  the  effects  of 
the  fuselage  on  the  antenna  performance  the  radiation  pattern  of  an 
electric  dipole  radiating  in  free  space  is  also  presented  in  Fig.  13. 
The  effect  of  the  aircraft  fuselage  on  the  radiation  patterns  is 
clearly  shown  through  the  comparison  between  these  two  patterns.  One 
should  expect  higher  level  of  sidelobes  in  the  forward  region  as  shown 
in  the  pattern  because  of  the  location  of  the  antenna  which  is  near 
the  front  portion  of  the  fuselage.  The  ripple  shown  in  the  pattern  is 
a  result  of  the  scattering  by  the  structure  of  the  fuselage  and  this 
is  taken  into  account  by  the  various  diffractions  among  the  wedges. 
Note  that  the  Imperfect  nulls  occur  at  0*90**  and  270**  should  be 
expected  because  of  the  un-symmetri cal  antenna  location  as  mentioned 
above.  Thus,  the  performance  of  the  antenna  is  greatly  affected  by 
the  aircraft  structure  (especially  the  fuseU,,.;  shape). 
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In  order  to  Improve  the  theoretical  simulation  of  an  aircraft 
fuselage,  a  curved  leading  edge  model  Is  Introduced.  This  model  Is 
shoMn  In  Fig.  14  and  has  a  hemisphere  replacing  the  cone  on  the  leading 
edge.  The  sphere  has  the  same  radius  as  that  of  cylinder.  This 
blunt  model  might  In  some  cases  be  a  better  approximation  of  an  actual 
fuselage  than  the  previous  two  models.  Creeping  wave  theory  Is 
employed  with  the  wedge  diffraction  technique  to  solve  this  problem. 

The  diffracted  fields  In  ray  form  from  edges  4,  5  and  6  are  the  same  as 
those  in  the  previous  case.  But  for  the  leading  edge,  the  creeping 
wave  expression  Is  now  applied  to  describe  the  lossy  surface  wave 
propagating  around  the  curved  surface  so  that  the  effect  of  the  scatter 
Ing  by  the  blunt  structure  could  be  properly  examined.  The  general 
form  for  a  diffracted  ray  radiating  from  point  Q2,  some  point  on  a 
curved  surface.  In  the  scattered  field  direction  with  the  source 
located  off  the  curved  surface  Is  given  by 


(71) 

where 


“cwWj)  =  W 


D„(Q,)DJQj,)  e 


-/a  (5-)dS^ 


■  the  point  on  the  curved  surface  at  which  the  direct 
ray  from  the  antenna  Is  grazingly  Incident. 

Qg  ■  arbitrary  diffracting  point  on  the  curved  surface. 

«  the  distance  or  arc  length  which  the  surface  ray 
travels. 

Rf(Q.|)  ■  the  Incident  ray  at  point  Q^. 

The  Incident  wave  excites  an  Infinity  of  surface  wave  inodes, 
as  the  summation  In  Eq.  (71)  Indicated,  each  with  Its  own  diffraction 
coefficient  and  attenuation  constant;  only  the  dominant  mode  Is  being 
considered  In  this  analysis  for  the  sake  of  simplicity,  I.e.,  m>0. 
Hence  Eq.  (71)  Is  reduced  to 
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-O-S,  -jks, 

(72)  R^iQ^)  »  Ri(Qi)Do(Qi)Dp(Q2)  e  °  ’  e  ^  . 

This  is  a  simplified  expression  for  a  diffracted  surface  ray  at  any 
arbitrary  direction  with  the  source  located  off  the  curved  surface. 

Note  that  the  elimination  of  the  integration  in  Eq.  (71)  is  due  to 
the  fact  the  attenuation  constant  is  not  a  function  of  ti ,  arc 
length,  but  a  constant  for  a  circular  cylinder  or  a  sphere  given  by[10] 


where  the  hard  surface  boundary  condition  is  employed  and  "a"  is  the 
radius  of  the  sphere.  The  diffraction  coefficient  D  is  given  as[10] 


(74) 


D  2  ^  ^  “■  (ka)l/3 

°h  A 


where  d^  is  equal  to  1.083. [10] 

Although  the  mounting  of  a  hemisphere  to  a  circular  cylinder 
produces  a  continuous  radius  of  curvature  and  first  derivative,  the 
second  derivative  of  the  radius  of  curvature,  however,  is  by  no  means 
continuous.  It  is  discontinuous  at  the  sphere-cylinder  junction.  Even 
so,  the  diffraction  at  the  discontinuity  Q-j  is  so  small  that  one  can 
neglect  its  effect  since  the  dimensions  of  the  model  are  quite  large 
in  terms  of  wavelength.  This  assumes  that  the  incident  ray  is  equal 
to  the  surface  ray  at  Qi.  Based  on  this  approximation,  the  diffraction 
coefficient  at  is  found  by 

(75)  R^(Q,)  -  F(e-0‘*)  e‘'^'“* 

(76)  R^^(Q^)  -  R^(Q,)  D^2(q^)  .  p(e.o«)  .  D^^CQ^)  . 

By  equating  these  two  equations,  D^(Q^)  is  given  as 
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.■9> 


(77)  D^{Ql)=1  • 

Therefore,  the  creeping  wave  expression  In  any  direction  can  be 
written  as 

(78)  R^(0)  F{  •O'*) 

-e[i^ka)^^^  (%)^^^  +  jka] 

e 

180"<e<360“ 

With  the  appropriate  phase  terms  Included  and  the  various 


regions  considered,  the  elevation  plane  radiation  pattern  Is  given 
by 

(79) 

*  R5(e)+Rg(e) 

O^0<gg 

(80) 

•  R5(e)+R5(e)+R5(e) 

65<e<180*-65 

(81) 

Rt(8)  •  RsC8)-*-R6(»)+R5(®)+’^4(®) 

18O'»-B5<0<18O‘* 

(82) 

RT(e)  • 

180*£«i»  SO^+gg 

(83) 

R|-(e)  -  Rg^(e)+R5(e)+R^(e) 

18O'+Bg<0;^36O*-B5 

(84) 

R|-(e)  -  R^(e)+R^(e) 

360*-B5<e<360* 

The  calculated  elevation  plane  radiation  patterns  for  a 
circumferential  slot  and  monopole  mounted  on  this  particular  model 
are  Illustrated  In  Figs.  15  and  16.  Those  with  cone  type  leading 
edges  are  also  presented  In  Figs  15  and  16.  From  these  figures  one 
observes  that  the  ripple  level  Is  decreased  considerably  for  the 


“jfjr 

s-  (ka)l/3, 


38 


blunt  model.  Thus  one  concludes  that  the  significant  effect  of  the 
curved  surface  Is  reducing  the  ripple  effect  and  smoothing  the  pattern. 

Summary  and  conclusion 

The  radiation  patterns  In  the  elevation  plane  for  an  antenna 
mounted  on  a  simple  fuselage  model »  composed  of  cones,  spheres,  and 
cylinders,  can  be  obtained  by  treating  this  three  dimensional  model 
as  a  tMo  dimensional  problem  through  the  use  of  modem  diffraction 
techniques.  The  results  obtained,  using  these  techniques  and  approx1> 
nations  compare  very  favorably  with  the  measured  results.  The 
structure  of  the  aircraft  fuselage  has  a  very  significant  effect  on 
the  antenna  performance  with  respect  to  distortion  of  antenna  patterns. 
For  exanple,  the  effect  of  a  cone  type  structure  often  Increases  the 
ripple  on  the  radiation  patterns  whereas  that  of  a  curved  surface 
tends  to  smooth  the  patterns.  The  location  of  the  antenna  Is  also 
an  Important  factor  In  determining  the  radiation  patterns  In  the 
elevation  plane.  In  addition,  the  pattern  In  the  Illuminated  region 
of  an  antenna  Is  simply  Its  direct  radiation  while  that  In  the  shadow 
region  Is  completely  dominated  by  the  scattering  of  the  structure 
of  aircraft  fuselage.  Although  these  simple  models  Illustrate  the 
effects  of  fuselage  structure  on  radiation  patterns,  they  fail  to 
Include  the  general  variety  of  aircraft  fuselages  encountered  In 
practice.  Thus,  a  more  general  representation  of  the  aircraft 
fuselage  Is  presented  In  the  next  section  in  which  any  convex  body 
(fuselage)  can  be  modelled  to  a  good  approximation. 

IV.  ELEVATION  PATTERNS  FOR  GENERAL  FUSELAGE  MODELS 

The  radiation  pattern  in  the  elevation  plane  for  an  aircraft 
fuselage  was  Initially  studied  using  a  basic  model  composed  of 
cylinders,  spheres,  and  cones  as  discussed  In  the  last  section.  The 
solutions  were  based  on  the  two-dimensional  analysis  which,  as 
demonstrated  by  Ryan, [9]  also  Is  valid  for  the  three  dimensional 
elevation  plane  patterns.  Some  of  the  results  of  this  stu4y  were 
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illustrated  in  Figs.  10  through  16.  Even  though  those  solutions  are 
valid  for  the  simple  models,  they  really  are  not  general  enough  to 
include  the  wide  variety  of  aircraft  fuselages  encountered  in  practice. 
For  example,  an  aircraft  fuselage  is  usually  a  convex  body  that  can 
not  be  completely  described  by  simple  mathematical  equations.  Hence, 
in  practice,  an  aircraft  is  often  specified  by  points.  Consequently 
a  new  approach  is  developed  in  which  a  set  of  discrete  points  is  used 
to  describe  the  shape  of  the  fuselage.  In  this  way  any  convex  fuse¬ 
lage  shape  can  be  included  in  our  general  solution  in  a  rather  simple 
manner. 

In  addition,  more  general  ellipsoidal  models  are  developed 
using  a  more  rigorous  "6TD"  solution.  This  “GTD"  solution  is 
developed  specifically  for  a  source  located  on  the  curved  surface, 
which  is  in  contrast  to  the  case  discussed  in  section  II  where  the 
source  was  located  off  the  curved  surface.  The  most  general  model 
being  considered  is  a  composite  (back  to  back)  ellipsoid,  however, 
elliptical  cross-section  cylinders  are  also  good  models  to  simulate 
a  fuselage.  In  order  to  illustrate  the  validity  of  this  new  method, 
called  "section  matching  method”,  the  radiation  patterns  for  circular 
cylinders  for  slot  and  monopole  antennas  are  first  calculated. 

Basically  this  section  is  written  to  illustrate  the  validity 
of  the  "Section  Matching  Method”.  The  only  geometries  considered  in 
its  application  are  for  those  that  can  be  solved  using  the  more 
rigorous  "GTD"  solutions.  In  the  future  this  method  will  be  used 
beyond  the  realm  of  the  "GTD”  solution  taking  into  account  more 
complex  fuselage  shapes.  Not  only  is  this  important  in  terms  of  the 
profile  of  the  fuselage  in  the  elevation  plane  but  also  in  terms  of 
the  cross-sectional  shape  in  the  roll  plane. 


A.  GTD  Solution  For  General  Fuselage  Models 

1 .  Circular  cylinder 

The  far  field  for  various  sources  such  as  slot  and  nonopole 
antennas  nounted  on  an  infinitely  long  perfectly  conducting  circular 
cylinder  has  been  analyzed  by  Kouyouniuian.[17]  In  this  solution  the 
reciprocity  theorem  is  employed  so  that  the  radiation  pattern  is 
derived  from  the  plane  wave  scattering  by  a  cylinder.  The  physical 
optics  approximation  for  the  surface  current  is  used  to  describe  the 
field  in  the  region  in  which  the  source  is  directly  illuminated  by 
the  incident  plane  wave  (the  lit  region).  The  Fock  current  approxi- 
mationClS]  is  used  for  the  penunbra  (transition)  region,  whereas,  in 
the  deep  shadow  region  of  the  cylinder,  the  geometrical  theory  of 
diffraction  solution  is  applied  to  give  the  far  field  expression.  It 
is  noted  that  a  launch  coefficient,  which  relates  the  creeping  wave 
(GTD  field)  to  the  actual  surface  field  or  current,  is  Introduced  so 
that  the  GTD  solution  can  be  properly  employed.  This  launch  factor 
is  deduced  from  the  exact  solution  for  the  surface  current  on  the 
circular  cylinder;  and,  in  Its  asymptotic  form,  involves  the  well-known 
Airy  function. 

The  far-zone  radiation  pattern  for  a  thin  slot  mounted  on  an 
infinitely  long  perfectly  conducting  circular  cylinder,  equivalent 

A 

to  a  magnetic  line  source  ky,  with  radius  "a”  as  shown  in  Fig.  17  can 
be  obtained  by  employing  the  aforementioned  GTD  asymptotic  expressions 
for  the  various  regions  and  is  given  by  Reference  [17] 


(85) 


1)  Lit  Region:  Geometrical  Optics  Description 

I  cos(e^-e) 

Ry(p.0)  -20  ® 

2)  Transition  Region: 


(86) 


Ry(p.0) 


9 


g*(Ci)  e 


-jkaii*. 


+  g*(C2)  e 
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REGION  I  :  LIT  REGION 
REGION  H-  TRANSITION  REGION 
REGION  m:  DEEP  SHADOW  REGION 


Fig.  17.  Antenna  mounted  on  a  circular  cylinder:  GTD  solution. 
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where 


(87a) 


{87b) 


<l»1  ■  l®o  ■  ®l  " 
*2  -  r  -  l«o  - 


M 


•  6  -  e. 


4  -  » 


•ir/2<^0<^3Tr/2 


The  functions  g*(e^)  and  g*(C2)  are  the  conjugate  of  g(5^)  and  g(e2)* 
respectively.  The  function  g(c)  is  the  so  called  Fock's  function  for 
hard  boundary  condition  (au/an  -  0)  and  is  given  as 


The  magnitude  and  phase  of  Fock's  function  g(x)  are  given  in 
Reference  [18]. 


(89) 


3)  Deep  Shadow  Region:  Surface  Ray  Description 

Ryb/e)  .  (-2.)  . 

*  f  h  >  -[oll(«)]Sl  -JkSi 

I.  A1(-5,)e  •  'e  ' 


mwl 


h  2  -jks-l 
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where 


$1  -  ai|»^ 


and  i|ii  and  \i)2  are  defined  in  Eq.  (87a)  and  q^^  Is  the  root  of 
Ai’(-?r  ),  that  is,  Ai'(-?r)  *  0.  The  diffraction  coefficient 

W  ^  Vi 

[Dj|j(a)]‘  and  the  surface  propagation  constant  otJ|,(a)  are  given  below, 
respectively. [15] 


[ol'(a)^  -  -57J — 
"  2^/6  ^ 


-i— 

,1/3  ,  ^12 


kax-2/3 


r\  0.1. 

l-air  •  ® 


a  6  /ka\l/3^ ,  .  fka\-2/3  i%  .  O.K  ' 

<lm 


The  roots  q„  and  A1(-q„)  can  be  found  in  Reference  [22].  It 
is  noted  that  the  surface  ray  expression  in  Eq.  (89)  can  be  rewritten 
in  the  following  way 


(92)  Rj,(p.9)  -  (-2.) 

k  9  ■y(«)So 

♦  [l{;(«)r  MKJ  e 

> 

where  Y(a)  ■  +  Jk>  propagation  constant  of  the  surface  ray. 

The  calculated  radiation  patterns  using  Eqs.  (85)  to  (92)  in 
the  xz-plane  (elevation  plane)  of  Fig.  17  for  various  cylinders  have 
been  shown  to  agree  very  favorably  with  those  obtained  from  the  modal 
solution,  especially  for  ka^S.O.  These  formulations  also  predict  the 
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radiation  pattern  fairly  well  for  small  cylinders  with  radii  (town 
to  the  range  from  1/3  to  1/4  wavelength. [19] 

2.  General  cross-section  cylintter 

The  results  for  the  circular  cylinder  are  extended  to  cylinders 
with  general  cross-sectional  shape  in  which  the  curvature  varies  along 
the  ray  path.  This  variable  curvature  has  a  strong  effect  on  the 
attenuation  constant  which  in  turn  affects  the  energy  propagation 
quite  significantly.  Thus,  the  representations  for  fields  in  the 
lit,  transition  and  deep-shadow  region  are  modified  to  include  the 
variable  curvature  effect  and  again  are  given  by 


1)  Lit  Region: 


(93) 


jkr,  cos(e-e^) 
Ry(e)  =  2  e  ^  0 


(94) 


2)  Transition  Region: 

P  P 

-j/  kds  -jr 

\9*Up)  e  ^  F(A)  +  g*(5g)  e 


Ry(e) 


III  I 

I 


kds 


F(B) 


3)  Deep-Shadow  Region; 


(95) 


#  a  * 


^(A)  AIKJ  . 


T^s)  ds 


’m' 


F(A)  +  dJ,(P)  dJ,(B)  Al(-q^)  e 


rJU)  d. 


F(B) 


where  A  and  6  are  the  points  that  the  surface  rays  launch  in  the 
direction  of  the  observation  point,  e;  r^  is  the  distance  from  the 
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origin  to  the  source  location;  and  Pg  is  the  radius  of  curvature  of 
the  curved  surface; 


and 


The  increment  of  the  arc  length  along  the  ray  trajectory  is  ds;  F(A), 

F(B)  are  the  phase  correction  factors  at  the  launching  points  A  and 
B;  g*(CB)  are  the  conjugates  of  and  gCSg).  Fock's  function 

as  defined  previously;  and  and  are  the  diffraction 

coefficients  at  the  source  location,  launching  points  A  and  B, 
respectively.  The  propagation  constant  is  y|J(s)  =  where 

is  the  ray  trajectory  attenuation  factor  for  the  hard  boundary  case 
and  k  Is  2ir/x.  In  addition,  AiC-q^^^)  is  the  Airy  function  evaluated  at 
-qV  where  ?  is  the  root  of  Ai'(-q„)  *  0.  Note  that  the  phase  reference 
point  for  these  formulations  is  at  the  origin.  Some  calculated  radiation 
patterns  in  xz-plbne  (the  elevation  plane)  enploying  these  formulations 
for  an  elliptical  cylinder  are  also  given  in  Reference  [19].  The  result 
compares  fairly  welt  with  the  exact  solution  obtained  by  Sinclair[20]  as 
shown  in  Fig.  18. 


The  same  technique  is  also  enployed  to  derive  the  general izea 
expressions  for  slot  antenna  with  soft  boundary  condition  and  in¬ 
finitesimal  monopole  antennjt.  The  equations  again  contain  the  Fock 
functions.  A  sunmary  of  the  far  field  expressions  for  the  circumfer¬ 
ential  slot,  axial  slot  and  Infinitesimal  monopole  are  presented  in 
Table  I. [17] 
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NORMALIZED  MAGNETIC  FIELD 


EXACT  (  SINCLAIR ) 
GTD 


•  •  ,  • 


SEMI  -  MAJOR  AXIS  :  a  >0.637  X 
SEMI>MINOR  axis:  b-0.390X 


160  180 


^  (DEGREES) 


Fig.  18.  Elevation  plane  patterns  with  a  circumferential 
slot  moulted  on  an  elliptical  cylinder. 

It  is  noted  that  the  oropagation  factor  in  the  far  field 
expression  has  the  form  for  the  circular  cylinder  case  and 

g~fr(s)ds  yjg  general  convex  surface  case.  Since  the  general 
formulations  can  not  be  calculated  in  a  closed  form  due  to  the 
complicated  expressions  y(s)  and/or  S,  the  results,  then,  are  often 
numerically  performed  via  a  modem  high-speed  digital  computer.  Hence, 
the  integration  is  often  performed  through  a  discrete  summation  process 
to  obtain  the  numerical  result,  i.e.. 


Y(a)  S  ^  ^  Y(a)  •  ASa  ■  y(«)  I  aS. 

i-1  ’  i-1  ' 


/y(s)  ds  ~  y(S^)  • 
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Based  on  the  fact  that  the  numerical  computation  In  a  digital 
computer  Is  a  discrete  process  and  that  an  aircraft  Is  generally 
specified  by  points  In  which  no  simple  equation  Is  available  to  give 
the  complete  description  of  the  aircraft  shape  as  mentioned  previously, 
a  new  approach  called  "the  section-matching  method"  has  been  developed. 
This  method,  by  definition.  Is  to  match  portions  of  a  general  surface 
through  the  use  of  a  polynomial  interpolation  process.  This  method 
has  proved  very  satisfactory  In  solving  for  the  radiation  pattern  for 
an  antenna  mounted  on  any  general  convex  body  such  as  an  aircraft 
fuselage  which  Is  our  main  Interest  In  this  report. 


B.  Section  Matching  Method 

In  order  to  predict  the  radiation  pattern  of  an  aircraft 
antenna  mounted  on  a  general  convex  fuselage,  the  section  matching 
method  can  be  applied.  In  this  method  the  continuously  curved  surface 
of  the  fuselage  Is  described  by  a  set  of  discrete  points.  Thus,  the 
generalized  6T0  formulations  as  shown  In  Table  I  now  have  to  be 
modified  so  that  they  can  be  used  for  the  discrete  points  case.  For 
example,  the  equations  for  a  circumferential  slot  are  easily  modified 
as  follows: 


(98) 


(99) 


1)  Lit  Region: 

jkr-cos(e-e_) 

2)  Transition  Region: 

A  -JkAS.  B  -jkAS 

Rje)  •  g*(LJ  F(A)  TT  e  ^  +  g*(eB)  •  F(B)  •  TT  e 
^  ^  1-1  ^  1-1 
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3)  Deep-Shadow  Region: 


(100) 


i-A 


+  dJ;(P)  D^B)  Ai{ 
m  m 


•  F(B)  • 


P  -jY"(S.)AS.l 
TT  e 
1-B 


where  A,  B  and  P  are  the  launching  points  and  source  point,  respectively; 
i  is  the  index  of  the  point  lying  between  the  launching  points  A  or  B 
and  source  point  P;  aS^  is  the  increment  of  ray  path  from  point  i  to 
point  (i-t-l):  dp  is  the  angle  in  the  direction  of  the  source  point 
refered  to  the  origin;  and  equals  [o|ji(S^)  +  jk]  and  is  the  propa¬ 

gation  constant  at  the  point  i  for  the  hard  boundary  condition. 

The  generalized  diffraction  coefficient  and  attenuation 
constant  as  presented  in  Reference  [15]  are  given  in  Table  II 
below.  Note  that  the  diffraction  coefficient  and  attenuation  con¬ 
stant  are  functions  of  the  longitudinal  radius  of  curvature  (p^) 
along  the  ray  trajectory,  its  first  and  second  derivatives  (p^  and 
with  respect  to  arc  length  along  the  ray  trajectory*  and  transverse 
radius  of  curvature  (p^p)>  In  this  two-dimensional  analysis,  p^^  Is 
assumed  to  be  infinite,  hence  the  terms  containing  ^ — In  the  expressions 
for  0.  and  a.  can  be  omitted.  It  is  obvious  that  the  logitudlnal 
radius  of  curvature  can  not  be  obtained  analytically  from  this  discrete 
set  of  points,  since  no  equation  is  given.  One  way  to  overcome  this 
difficulty  is  to  match  sections  of  the  fuselage  surface  using  a  numerical 
interpolatpry  approximation  in  which  an  (n-l)th  order  polynomial, 
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which  passes  through  n  points  In  a  given  section.  Is  determined.  In 
other  words,  a  continuous  polynomial  function  Is  found  numerically  to 
approximate  the  actual  curve  so  that  at  the  points  x^,  1*1,  2,  •••,  n, 
the  approximation  function  agrees  with  the  actual  function  f(x).  Once 
this  polynomial  Is  obtained,  the  radius  of  curvature  can  be  determined 
without  difficulty.  The  first  and  second  derivatives  of  the  radius 
of  curvature  can  also  be  obtained  simply  by  numerically  differentiating 
with  respect  to  the  arc  length.  Hence,  the  diffraction  coefficient 
and  attenuation  constant  are  known  and  Eqs.  (98),  (99),  and  (100)  can 
be  employed  to  calculate  the  radiation  patterns. 

In  order  to  determine  the  required  polynomial  y(x),  the  exact 
Lagranglan  Interpolation  formu1a[20]  Is  used.  To  Interpolate  n 
discrete  points,  the  Lagranglan  polynomial  y(x)  should  be  of  degree 
less  than  n  with  the  prcnerty  y(a^)  •  f(a^),  1*0,  ••*,  n-1,  where  f(x) 
Is  the  original  function  or  ctual  curve  and  a^  are  the  discrete 
points.  Thus,  the  polynomial  y(x/  can  be  written  as 

(101)  y(x)  *  Cg+C|X+C2X^+***+c^_^x'''"^  . 

At  the  discrete  points  a^,  1*0,  •••,  n-1, 

(102)  y(a^)  *  *  ^^*1^ 

or 

(103)  f (a^)  *  Vc,ap+C2ao^+C3a^^- •  •+c„.,ap"‘’ 
f(a,)  *  VCia,+C2a,^C3a,^...+c„.,a,"*’ 
f(e2)  ■  ®o^®l*2'’^2*2^''‘®3*2^'*“ ‘“'^0-1*2"  ^ 


where  f(af)  Is  known. 
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This  linear  system  of  equations  can  be  solved  by  talcing  the  Inverse  of 
matrix  A.  In  which  case 

(107)  c  *  A“^y 

Once  the  unknowns  c^  are  found*  the  polynomial  y(x)  Is  determined 
as  In  Eq.  (101).  The  radius  of  curvature  may  be  obtained  from  these 
equations  us1ng[21] 


The  ambiguous  sign  In  Eqs.  (108)  and  (109)  Is  chosen  so  that 
Pg  Is  always  positive  since  the  curve  is  convex. 
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Further,  the  first  and  second  derivatives  of  the  radius  of  curvature 
as  required  in  the  expressions  of  diffraction  coefficient  and  at¬ 
tenuation  constant  are  calculated  by  nunerically  differentiating  the 
radius  of  curvature  with  respect  to  the  arc  length  s.  Hence 

Eqs.  (98)  to  (100)  can  now  be  applied.  Consequently,  the  radiation 
patterns  for  an  antenna  mounted  on  a  general  convex  fuselage  can  be 
obtained. 

In  order  to  illustrate  the  accuracy  of  this  solution,  the 
radiation  patterns  for  an  antenna  mounted  on  circular,  elliptical, 
and  composite  elliptical  cylinders  are  calculated  using  this  method. 
First  these  cylinders  are  described  by  discrete  points  which  are 
generated  by  equations  for  a  circle  and  ellipse.  The  points  are 
generated  in  such  a  way  that  the  angle  difference  between  the  tangent 
direction  at  any  two  adjacent  points  is  less  than  2  degrees.  This 
is  due  to  the  fact  that  the  surface  ray  sheds  energy  in  the  tangent 
direction  and  a  reliable  radiation  pattern  is  likely  to  be  obtained  only 
If  computed  at  intervals  less  than  or  equal  to  2  degrees.  The  shedding 
angle  (tangent  angle)  at  any  point  on  the  fuselage  can  be  found 
approximately  by  the  following  relation  (see  Fig.  19). 

(110)  0  »  o  -  0. 

where 


o  »  tan 


h^i  -  yjl 

^i+i  - 


2pgi 

and  p^^  is  the  radius  of  curvature  at  point  1.  However,  p^^  is  an 
unknown  quantity  but  it  can  be  found  easily  as  will  be  discussed  later. 
Therefore,  one  should  use  a  trial  and  error  method  to  generate  and 
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Fig.  19.  Determination  of  shedding  angle. 

test  the  discrete  set  of  points  from  the  analytic  equations  In  order 
to  obtain  a  sufficient  number  of  points  to  describe  the  model.  After 
obtaining  the  required  discrete  points,  the  radius  of  curvature 
and  its  first  two  derivatives  and  p^  are  found  according  to  the 
method  just  discussed. 

A  Lagranglan  Interpolation  polynomial  of  degree  n-1,  where  n 
Is  chosen  to  be  an  odd  number  for  our  application.  Is  used  to  ap¬ 
proximate  the  function  which  passes  thr  .^h  any  n  sequential  points 
on  the  surface.  In  other  words,  one  Is  trying  to  match  a  section, 
which  contains  these  n  points,  of  the  geometry  of  particular  Interest 
(fuselage)  by  an  (n-l)th  order  polynomial  approximation.  After  the 
desired  polynomial  is  obtained,  Eq.  (108)  or  Eq.  (109)  Is  eeployed  to 
find  the  radius  of  curvature  at  the  center  point  of  this  section.  Tire 
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same  procedure  1s  applied  and  continued  until  the  radius  of  curvature 
at  every  point  on  the  surface  Is  found.  It  Is  noted  that  the  degree  of 
Lagranglan  polynomial  must  be  n-1  or  less  for  this  Interpolation  to  be 
exact,  I.e.,  the  error  at  these  n  points  Is  zero.  Note  also  that  the 
computer  time  required  to  solve  n  linear  equations  Is  proportional 

3 

to  n  .  Using  the  Crout  method  to  solve  the  n  x  n  matrix,  the  computer 

2 

time  Is  reduced  and  proportional  to  n  .  It  takes  approximately  one- 
fourth  the  tine  to  solve  n/2  equations  than  It  does  to  solve  n 
equations.  Thus,  for  the  sake  of  computer  time.  It  Is  wise  to  choose 
the  number  n  as  small  as  possible.  However,  the  accuracy  with  which 
the  polynomial  function  approximates  the  actual  surface  also  should 
be  taken  Into  careful  consideration.  Since  the  geometry  of  our 
Interest  Is  quite  smooth,  high-order  polynomials  are  not  necessarily 
required.  A  polynomial  of  degree  4  (or  n=5)  Is  thus  chosen  to  Inter¬ 
polate  th?  sections  of  our  present  fuselage  models. 

To  Illustrate  the  accuracy  of  our  section  matching  Interpolation, 
the  calculated  results  of  radius  of  curvature  and  Its  first  two  de¬ 
rivatives  from  a  set  of  discrete  points  generated  from  a  composite 
ellipse  with  two  different  senri-elllpses  placed  back  to  back  Is 
presented  In  Table  III.  Due  to  the  symmetry,  only  half  of  the  data 
is  shown.  Notice  that  the  radius  of  curvature  calculated  from 
discrete  points  compares  well  with  p  of  the  analytic  equation  of  the 
ellipses  except  at  0=90®  where  the  two  ellipses  are  joined  together. 

The  Inaccurate  data  for  the  derivatives  around  e*90®  and  e«270®  resulted 
from  the  error  accumulation  due  to  the  discontinuous  radii  of  curvature 
used  In  the  differentiations.  In  practice,  this  data  can  be  set  to 
equal  to  zeros,  since  Its  effect  Is  so  small  that  one  can  Ignore  It. 

Using  these  results  for  the  radius  of  curvature  and  its  derivatives 
along  with  Eqs.  (98),  (99),  and  (100),  the  radiation  pattern  in  the 
elevation  plane  for  various  antennas  mounted  on  circular  cylinders  with 
radius  a=lx  and  a«5x  are  calculated  and  presented  In  Figs.  20  to  25. 

Only  the  dominant  creeping  wave  terms  with  the  first  two  modes  are 
considered  In  these  patterns.  Notice  that  excellent  agreement  Is 
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TABLE  III 

COMPOSITE  ELLIPSE  (a=2.0  b=0.5  c»4.0) 


0“ 

Pc 

Pd 

i  • 

Pc 

Jd 

p'c 

Pd 

i 

0 

0.125000 

0.125000 

0.0 

0.00 

22.5 

22.493027 

5 

0.535155 

0.535156 

-3.507460 

-3.506249 

10.832942 

10.840757 

10 

1.829328 

1 .829360 

-5.298786 

-5.296271 

3.087004 

3.099108 

15 

3.385806 

3.385434 

-5.611507 

-5.602857 

-0.518659 

-0.453016 

20 

4.680069 

4.680116 

-5.250263 

-5.248033 

-2.413772 

-2.402165 

25 

5.624736 

5.6247^9 

-4.684848 

-4.680745 

-3.501232 

-3.535900 

30 

6.291131 

6.2912UI 

-4.102222 

-4.100489 

-4.169641 

-4.178182 

40 

7.100539 

7.100552 

-3.078533 

-3.077841 

-4.898717 

-4.893385 

50 

7.530519 

7.530523 

-2.260492 

-2.260237 

-5.255688 

-5.253419 

60 

7.771503 

7.771510 

-1 .590659 

-1 .590547 

-5.447772 

-5.446614 

70 

7.907789 

7.907797 

-1.015261 

-1 .015208 

-5.554036 

-5.553343 

80 

7.978184 

7.978191 

-0.494959 

-0.494937 

-5.608276 

-5.607783 

90 

7.999997 

12.800099 

-0.000001 

-824.250901 

-5.624997 

11834.8320 

100 

31.977036 

31.977060 

0.520462 

0.520457 

-5.901922 

-5.901795 

no 

31 .902405 

31 .902448 

1 .072628 

1 .072616 

-5.887842 

-5.887655 

120 

31.755432 

31.755494 

1.6%152 

1 .696123 

-5.860015 

-5.85%94 

130 

31.487198 

31 .487221 

2.451000 

2.450929 

-5.808909 

-5.808264 

140 

30.979645 

30.979704 

3.442989 

3.442778 

-5.711161 

-5.709542 

150 

29.907730 

29.907842 

4.885924 

4.885040 

-5.499908 

-5.494330 

160 

27.146667 

27.147224 

7.257616 

7.250326 

-4.922618 

-4.900683 

170 

17.579269 

17.581758 

11.146346 

11.122554 

-2.387525 

-2.318478 

175 

6.324325 

6.301114 

11.098542 

10.779890 

3.471679 

2.321958 

180 

0.062500 

0.062500 

0.000119 

47.250000 

47.241440 

60 


RELATIVE  POWER  (48) 


Fig.  21.  Elevation  plane  patterns  with  an  axial  slot 
mounted  on  a  circular  cylinder. 
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Fig.  22.  Elevation  plane  patterns  with  aiLlnflnitesIraal 
monopole  mounted  on  a  circular  cylindfer. 


Fig.  23.  Elevation  plane 
slot  mounted  on 


RELATIVE  POWER  (dB)' 


Fig.  25.  Elevation  plane  patterns  with  an  infinitesimal 
monopole  mounted  on  a  circular  cylinder. 
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observed  between  these  patterns  and  those  obtained  using  the  modal 
solution.  These  results  Illustrate  the  validity  of  the  section  matching 
method  for  this  simple  model. 

This  section  matching  method  now  Is  extended  to  obtain  radiation 
patterns  In  the  elevation  plane  for  a  simulated  fuselage  model  with  an 
elliptical  profile.  The  calculated  radiation  patterns  for  general 
antennas  such  as  circumferential  slot,  axial  slot  and  Infinitesimal 
monopole  mounted  on  an  Infinitely  long  elliptical  cylinder  with  a 
setiri -major  axis  a=0.637x  and  semi -ml  nor  axis  b=0.390x  are  shown  In  Figs. 
26  to  28.  The  comparison  between  the  continuous  GTD  solution  and  the 
section  matching  method,  again.  Is  very  satisfactory.  In  addition 
this  latter  method  has  been  applied  to  more  general  fuselage  models 
such  as  composite  elliptical  cylinders.  The  radiation  patterns  In 
the  elevation  plane  for  a  circumferential  slot  mounted  at  various 
locations  on  a  composite  elliptical  cylinder  are  shown  In  Fig.  29.  The 
results,  again,  compare  very  favorably  with  the  more  rigorous  GTD 
solution  [19].  The  elevation  plane  patterns  for  axial  slot  and 
Infinitesimal  monopole  are  also  shown  In  Figs.  30  and  31,  respectively. 
Based  on  these  results,  the  pattern  In  the  Illumination  region  Is 
simply  the  direct  radiation  by  the  antenna  Itself.  On  the  other  hand, 
the  pattern  In  the  shadow  region  Is  completely  dominated  by  the 
fuselage  structure  through  the  diffraction  mechanisms.  Thus,  the 
complete  radiation  pattern  for  an  antenna  mounted  on  a  fuselage  Is 
greatly  affected  by  the  fuselage  as  one  should  expect.  The  location 
at  which  the  antenna  Is  mounted  Is  also  an  Inportant  factor  In  pre¬ 
dicting  the  radiation  pattern.  Therefore,  one  should  foresee  a  better 
forward  coverage  for  slot  antennas  located  at  the  front  portion  of 
the  fuselage.  The  radiation  patterns  In  the  elevation  plane  fbr  a 
smaller  composite  elliptical  fuselage,  such  as  ai^l.Sx  and  a2‘‘2.0x 
and  b*0.5x,  are  presented  In  Fig.  32.  From  the  comparison  between 
Figs.  29a,  30a,  31a,  and  32  one  notices  that  there  are  fewer  backlobes 
but  more  energy  radiated  In  the  shadow  region  for  the  smaller  cylinder. 
The  larger  cylinder  has  more  backlobes  at  a  lower  average  energy 


67 


RELATIVE  POWER  (dB) 


Fig.  28.  Elevation  plane  pattern  with  an  Infinitesimal 
monopole  mounted  on  an  elliptical  cylinder. 


Fig.  29(a).  Elevation  plane  pattern  with  a  circumferential 
slot  mounted  on  a  conposite  ellipse  model. 
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Fig.  30(a).  Elevation  plane  pattern  with  an  axial  slot 
mounted  on  a  composite  ellipse  model. 
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Fig.  30(b).  Elevation  plane  pattern  with  an  axial  slot 
mounted  on  a  composite  ellipse  model. 
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Fig.  30(c).  Elevation  plane  pattern  with  an  axial  slot 
mounted  on  a  composite  ellipse  model. 


Fig.  32(a}.  Elevation  plane 
slot  moimted  on 


Fig.  32(b).  Elevation  plane  pattern  with  an  axial  slot 
mounted  on  a  composite  ellipse  model. 
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Fig.  32(c).  Elevation  plane  pattern  with  an  Infinitesimal 
monopole  mounted  on  a  composite  ellipse  mo^l 


level.  This  Is  simply  due  to  the  size  of  the  fuselage  and  the  inter¬ 
action  between  the  two  surface  waves  propagating  in  the  opposite 
directions  from  the  source  In  which  case  the  two  ray  paths  play  a 
most  Important  role.  That  Is,  the  longer  the  ray  path  is,  the  more 
the  surface  ray  attenuates.  Again,  this  shows  the  significant  effect 
which  the  structure  of  the  aircraft  fuselage  plays  in  terms  of  the 
antenna's  performance.  Therefore,  an  antenna  designer  should  take 
the  structure  of  aircraft  into  careful  consideration  in  the  design 
of  airborne  antennas  in  order  to  achieve  his  specific  needs  for  a  given 
application. 

V.  CONCLUSIONS 

Modem  diffraction  theory  in  the  form  of  the  geometrical  theory 
of  diffraction  (GTO)  and  its  extensions  including  wedge  diffraction 
theory  and  creeping  wave  theory  is  employed  to  analyze  the  elevation 
patterns  of  airborne  antennas  mounted  on  the  aircraft  fuselage.  Some 
simple  models  composed  of  cylinders  and  cones  are  first  used  to 
simulate  the  fuselage  so  that  wedge  diffraction  theory  could  be 
applied  to  calculate  the  radiation  patterns  in  the  elevation  plane 
for  antennas,  such  as  a  circumferential  slot,  an  axial  slot  and  a 
radial  monopole,  mounted  on  the  fuselage.  Based  on  this  study,  it 
was  found  that  the  fuselage  shape  has  a  dominant  effect  on  the  elevation 
plane  radiation  pattern.  The  radiation  patterns  calculated  through 
these  models  compare  very  favorably  with  the  measured  results.  In 
order  to  Introduce  the  curvature  of  an  actual  fuselage,  the  theoretical 
models  are  modified  to  Include  curved  surfaces.  A  combination  of 
wedge  diffraction  techniques  and  creeping  wave  theory,  then,  was  applied 
to  solve  the  new  geometry.  Again,  the  effect  of  fuselage  Is  clearly 
shown  In  the  radiation  patterns. 

In  order  to  Include  a  more  general  geometry  for  a  fuselage, 
a  new  approach  has  been  developed  In  which  the  fuselage  profile  Is 
described  by  a  set  of  points.  With  this  and  a  more  rigorous  "GTD" 
solution,  the  radiation  pattern  In  the  elevation  plane  fbr  an  antenna 


mounted  on  the  fuselage  can  be  calculated.  This  method,  called  “the 
section  matching  method,"  is  essentially  a  numerical  interpolation  in 
sections  described  by  a  set  of  points  which  defines  the  profile  of 
the  aircraft  fuselage.  With  this  interpolating  technique  the  radius 
of  curvature  and  its  first  two  derivatives  at  the  discrete  points 
are  found  nunerically.  Consequently,  the  generalized  “GTD"  solution 
which  is  a  function  of  the  radius  of  curvature  of  the  fuselage  can 
be  modified  to  obtain  the  radiation  patterns  for  a  fuselage  outlined 
by  discrete  points. 

In  order  to  illustrate  the  accuracy  of  this  solution,  the  section 
matching  method  is  applied  to  circular,  elliptical  and  composite 
elliptical  cylinders  in  which  the  composite  elliptical  cylinder  is 
considered  to  be  a  convenient  model  for  simulating  an  aircraft  fuselage. 
With  the  radiation  patterns  calculated  from  this  method  and  a  modal 
solution,  excellent  agreement  is  achieved  for  the  circular  cylinder 
case.  The  results  for  elliptical  and  composite  elliptical  cylinders 
also  compare  very  favorably  with  the  more  rigorous  continuous  GTD 
solutions.  These  good  results  verify  the  validity  of  this  section 
matching  method.  Based  on  this  study  one  concludes  that  the  pattern 
in  the  illuminated  region  is  simply  the  direct  radiation  by  the 
antenna  itself  whereas  that  in  the  shadow  region  is  completely 
dominated  by  the  fuselage  structure  through  the  diffraction  mechanisms. 
Thus,  the  aircraft  structure  has  a  significant  and,  often,  dominant 
effect  on  the  antenna  system  performance  especially  when  the  fuselage 
is  strongly  illuminated.  Therefore,  in  designing  antenna  systems  fbr 
airborne  applications,  the  antenna  designers  should  take  the  aircraft 
structure  into  careful  consideration  in  order  to  achieve  their 
specific  objectives. 
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The  radiation  patterns,  in  the  elevation  plane,  of  airborne 
antennas  mounted  on  the  aircraft  fuselage  have  been  analyzed  by 
modem  diffraction  theory  including  wedge  diffraction  theory  and 
creeping  wave  theory.  It  is  found  that  the  fuselage  shape  has  a 
dominant  effect  on  the  elevation  patterns.  In  order  to  analyze  a 
general  convex  shape,  a  new  approach,  "section  matching  method", 
has  been  developed  in  which  ^e  fuselage  profile  is  described  by  a 
set  of  points. 

This  new  method  is  applied  to  some  general  geometrical  shapes 
which  approximate  aircraft  fuselage.  The  composite  elliptical 
cylinder  is  the  most  general  model  considered.  The  validity  of  this 
method  has  been  verified  by  modal  and  more  rigorous  GTD  solutions. 
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